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SECTION  I 


INTRODUCTION 

Statistical  data  accumulated  from  five  nations  has  brought  to  light 
the  devastating  effects  of  high-speed  ejections  from  aircraft  (Glaister, 
1975).  That  is,  the  data  show  that  windblast  forces  increase  with  aircraft 
speed  to  the  point  where  an  otherwise  five  or  ten  percent  limb-flail  injury 
rate  rises  to  40  percent  or  more.  This  is  far  from  negligible,  but  control 
of  the  forces  that  produce  excessive  motion  of  the  limbs  of  an  ejection 
seat  occupant  can  only  be  achieved  if  we  increase  our  ui  derstanding  of  the 
aerodynamic  loading  to  which  a pilot  is  exposed  during  high-speed  ejections. 
In  a recent  report  (Schneck,  1976)  it  was  shown  that  it  is  feasible  to 
formulate  a mathematical  model  which  can  predict  such  aerodynamic  loading. 
Results  from  this  model  can  then  be  incorporated  into  the  Aerospace  Medical 
Research  Laboratory's  modified  Calspan  Model  of  the  Articulated  Total  Body 
(ATB)  in  order  to  assess  the  kinematics  of  limb  motion  under  the  action  of 
specified  aerodynamic  forces.  In  the  earlier  report,  potential  flow  solu- 
tions were  presented  for  estimating  the  pressure  distribution  around  the 
forearm  of  a human  body  subjected  to  windblast.  Work  completed  under  the 
present  grant  has  expanded  these  preliminary  results,  and  extended  them  in 
order  that  future  studies  may  include  some  effects  of  flow  separation. 

The  body  model  used  here  for  the  analysis  of  windblast  forces  is  repre- 
sented again  by  the  15-linkage  system  of  spherical,  circular-cylindrical, 
truncated-conical  and  flat-plate  segments  shown  in  Figure  1 (see  Schneck, 
1976,  for  a complete  description).  Moreover,  specific  attention  is  focussed 
on  the  events  which  take  place  where  a body  segment  is  either  in  contact 
with  a restraining  surface,  such  as  the  forearm  in  contact  with  an  arm  rest, 
or,  equivalently,  where  two  body  segments  are  in  contact  with  ona^anether, 

ri 


such  as  the  upper  arm  pressing  against  the  thorax.  For  this  situation,  the 
complex  velocity  potential  for  the  aerodynamic  cross-flow  over  the  limbs  (in 
the  absence  of  flow  separation)  has  been  shown  to  be  given  by: 

3.7T 

w = TraV  coth  — , (1) 

X 

which  corresponds  to  the  cross-flow  streamline  pattern  depicted  in  Figure  2 
for  the  cross-sectional  geometric  configuration  shown,  with  x = y + iz  and 
w = 4>  + if.  The  streamlines  illustrated  are  non-dimensionalized  with  respect 
to  V and  a. 

In  the  section  which  follows,  the  forces  tending  to  dislodge  the  limbs 
from  one  another,  or  from  a restraining  surface,  are  calculated  in  the  absence 
of  flow  separation,  i.e.,  for  the  situation  depicted  in  Figure  2.  Section 
III  then  describes  how  the  results  can  be  modified  to  include  the  more 
realistic  physical  effects  of  separation  of  flow  around  the  blunt  body 
segments . 


SECTION  II 

LIMB  DISLODGING  FORCES  IN  THE  ABSENCE  OF  FLOW  SEPARATION 

Let  p designate  the  aerodynamic  pressure  existing  at  a point  along  the 
cylindrical  surfaces  illustrated  in  Figures  2 and  3.  If  Y defines  the 

Z~cL 

azimuthal  coordinate  such  that  tan  y = (non-dimensionalized  radius  a = 1 

y 

in  Figure  2),  then  the  differential  surface  element  a(dy)  having  unit  length 
along  the  axis  of  the  limb  has  a force  on  it  of  magnitude  df  = p(a)(dy) 
directed  radially  inward  (see  Figure  3) . The  component  of  this  force  tend- 
ing to  push  the  cylinders  together  (or  pull  them  apart)  is  thus  given  by 

df  * -(p  sin  y)(a  dy) , so  that  the  net  vertical  force  becomes: 
z 

3v/2 

f * J -ap  sin  y dy  per  unit  length 
Z -ir/2 


(2) 


Now,  the  pressure,  p,  was  obtained  earlier  (Schneck,  1976)  in  the  form  of  a 
pressure  coefficient  which,  for  the  cylindrical  configuration  illustrated  in 
Figure  2,  is  given  by: 


Before  substituting  equation  (3)  into  equation  (2)  and  performing  the  indi- 
cated integration,  it  is  convenient  to  make  some  transformations  and  redefine 

certain  variables.  Referring  to  Figure  3,  we  note  that  the  equation  of  the 
2 

circle  is  r = 2az,  while  r = 2a  sin  0 and  y = r cos  0.  Thus, 


2 . 2 . 2 . „ 

Furthermore,  sin  y = — — = - — 1 = — — 1 = a-s  3 1 = 2sin  0-1 

a 3 2a2  2a2 

2 2 2 2 
= sin  0 + (sin  0 - 1)  = sin  0 - cos  0 = - cos  20 


so,  y = 20  + — j-  , and  dy  * 2 d0. 


2 4 2 4 4 4 4 

an  air  air  ir 

Finally,  j * —5-  = 5 r—  5 2~ 

4z  4r  16a  sin  0 16  sin  0 sin  0 


ll 1 d6 

8 2 d 0 (6) 

sin  0 


2 2 2 
And,  sin  0 - cos  0 , 2 , 46 

2 = 1 - cot  0 = 1 2~ 

sin  0 ir 

Putting  equations  (3)  through  (7)  into  (2)  gives:  (Limits  are  on  6) 


— 00  M 

/-(a)  »1  ,2  2 -,r3  1 d*  . 4 . 

jy  pUo  Sln  0 1 " ( ~B>  TTT  d^  SGCh  6 

\ sin  0 


||  sech4  6 + pc  j [-  cos  20] 2 d0  (8) 


o 2 2 12  2 

Now,  /(-  f pUQ  sin  a + pQ][-cos  20 ] 2 d0  = 0,  so,  with  qQ  = y pUQ  sin  a, 

-ir 
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and  the  use  of  equation  (7),  equation  (8)  simplifies  to: 

3 j u/.2  +<=°  , 2 

f ~ f 2(aq  sech  6 [ — — — 1 ] d 6 = -aq  sech  6 [6  - ^7-]  d6  (S 

z , o o 2 o 4 

+00  Tf  -CD 

Equation  (9)  contains  integrals  which  are  special  cases  of  the  general  form 
(see  Gradshteyn  and  Ryzhik,  page  124) : 


m , „ m-1  , , . m , , 

6 d6  _ m6  cosh  6 + (n-2) 6 sinh  6 


cosh  6 


(n-1)  (n-2) cosh11  1 6 

, , \ m-2 , _ m , 

_ m(m-l)  y 6 d<5  + n-2  y da 

(n-1)  (n-2)  cosh11  2 5 n-1  cosh11  2 6 


In  the  first  integral  of  equation  (9),  m = 2 and  n = 4.  In  the  second  inte- 
gral, n is  again  equal  to  4,  but  this  time  m = 0.  Noting  that  the  integrals 
are  symmetric,  we  may  write,  for  the  case  m = 2: 

2 CO  00  00  2 

/ 52  sech4  6 d 6 = 6-C-*S-h  6 + 6 sinh^^  _ 1 f _ji|_  +2  f _j_d^  (11) 

0 3cosh  6 o cosh  ,5  o cosh  5 

The  first  term  on  the  right  hand  side  of  equation  (11)  evaluates  to  zero  by 

00 

application  of  L'Hospital's  rule.  The  second  term  integrates  to  tanh  6]^  = 1 
(see  Gradshteyn  and  Ryzhik,  page  99)  and  the  third  integral  is  of  the  form 

00  (22k  _ 2->TT 2k 

(Gradshetyn  and  Ryzhik,  page  353):  f x = xr |b~,  |,  which, 

o coshZ  b6  b(2b) 

v2  1 

for  k=l=b  evaluates  to  yj  with  the  Bernoulli  number  | | “ ^ • Putting  all 

of  this  together,  the  first  integral  of  equation  (9)  may  be  evaluated  in 

2 

closed  form  to  yield  ( tt  - 6)/9.  Now,  for  the  case  m = 0,  equation  (10) 

becomes:  S sech4  6 d6  =»  8^nh^'5  + y J"  — (12) 

o 3cosh  6 0 cosh  6 

As  before,  the  first  term  on  the  right  hand  side  of  equation  (12)  evaluates 

to  zero,  and  the  second  term  integrates  to  unity.  Thus,  the  second  integral 

2 

of  equation  (9)  may  be  evaluated  in  closed  form  to  yield  (4/3) (-n  /4) . We 


may  then  write: 


10 


(13) 


- ,ti2-6  4tt  . aqo1T  , „ 2 , 

fz  = -aV  [~ T - "12  ] = ~ f-2lT  - 61  = 9aq0 

Note  the  positive  resultant  sign  of  equation  (13),  indicating  that  the  net 
force  per  unit  length  of  limb  is  in  the  positive  z-direction,  or  tending  to 
cause  dislodgment  of  the  limb  from  a restraining  surface  (which  could  also  be 
another  limb) . 

In  order  to  evaluate  equation  (13)  further,  we  need  to  say  something 

about  qQ  and  a.  As  a reasonable  approximation,  we  may  begin  by  assuming 

that  the  radius  of  the  human  forearm  is  on  the  order  of  1.5  inches  (average) 

and  that  the  length  of  this  portion  of  the  arm  is  about  one  foot.  Thus,  the 

total  force  acting  on  the  forearm  can  be  approximated  by  the  relation 
9 

f = -T-  q . Now,  consider  an  ejection  taking  place  at  an  altitude  of  10,000 
t o o 

feet.  At  this  height,  the  density  of  air  relative  to  a known  standard  is 

given  by  (Shames,  page  540)  = 0.7385,  where  p = 0.002378  slug/ft  , and 

o 

the  speed  of  sound,  c = 1,078  feet  per  second.  Thus, 


qo  " 2 


= i p [U  ^/c^]c^sin^  a = 1020M^sin^  a, 


p oo 
o 


where  M corresponds  to  the  Mach  Number  of  the  flow.  Substituting  this  result 
into  the  equation  for  f gives  the  total  limb-dislodging  force  as  a function 
of  Mach  Number  and  Angle  of  Attack  at  the  specified  altitude  of  10,000  feet. 


The  final  equation  is: 


f = 1148  sin^  a 


(14) 


Figure  4 shows  equation  (14)  plotted  as  limb  dislodging  force  vs.  Mach  Number 
with  the  angle  of  attack  appearing  as  a parameter.  Although  the  figure  illus- 
trates force  behavior  in  the  supersonic  range  as  well  as  the  subsonic,  it  is 
recognized  that  the  results  past  M = 1 are  somewhat  unrealistic  in  that  the 
analysis  has  thus  far  neglected  the  presence  of  any  shock  waves  in  the  flow. 

In  any  case,  the  subsonic  values  for  f are,  in  themselves,  very  revealing. 

For  example,  at  the  higher  angles  of  attack  one  observes  a rapidly  increasing 
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limb-dislodging  force.  The  significance  of  these  large  forces  may  be  appre- 
ciated by  examining  them  in  relation  to  the  super-imposed  ordinate  axis  which 
is  a measure  of  several  published  probability  results.  That  is,  in  a study 
conducted  by  Horner  and  Hawker  (1973)  it  was  found  that  a pilct's  average 
grip  retention  is  such  that  the  probability  of  his  letting  go  is  100%  if  the 

dislodging  force  exceeds  some  600  pounds.  From  Figure  4,  we  see  that  f = 600 

o o 

pounds  when  M = .72  for  a = 90  and  when  M = .83  for  a = 60  . Thus,  at 
these  higher  angles  of  attack  a pilot's  musculo-skeletal  system  is  not  likely 
to  withstand  the  tendency  for  dislodgment  from  a restraining  surface  if  he 
is  ejecting  at  Mach  Numbers  in  excess  of  around  0.7.  Similarly,  in  a study 
conducted  by  Payne  (1975)  it  was  found  that  the  probability  of  major  flail 
injury  is  around  100%  if  the  ejection  Mach  Number  exceeds  1.25.  This  would 
appear  to  be  self-consistent  with  the  predicted  limb  dislodging  forces  as 
interpreted  above  together  with  the  results  of  Horner  and  Hawker.  To  this 
level  of  approximation,  therefore,  one  would  have  to  conclude  that  the  genera- 
tion of  stagnation  points  in  the  flow  produces  forces  that  can  cause  limb- 
dislodgment  (with  subsequent  flail  and  possible  serious  injury).  Moreover, 
these  forces  are  sensitive  to  the  angle  at  which  the  limb  intercepts  the 
flow,  such  that  the  higher  the  angle,  the  greater  the  tendency  for  dislodgment. 
And  finally,  the  forces  increase  rapidly  with  speed  of  ejection,  which  corre- 
lates well  with  the  finding  that  windblast  injuries  increase  dramatically  as 
a function  of  airspeed  (Glaister,  1975). 

Figure  5,  taken  from  Payne,  Hawker  and  Euler  (1975)  is  a photograph  of 

o o 

an  ejection  seat  occupant  sitting  in  the  ACES-II  Seat  at  -15  Yaw  and  -15 
Pitch,  during  a wind-tunnel  simulation  of  an  ejection  from  the  F-105.  The 
arrows  point  out  several  critical  regions  of  the  pilot-ejection-seat  configura- 
tion where  flow  stagnation,  with  its  consequent  limb-dislodging  force  distribu- 
tion, is  likely  to  occur.  Observe,  in  particular,  those  regions  where  the 
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pilot  is  gripping  the  seat,  where  his  upper  arms  are  in  contact  with  his 
torso,  where  the  lower  legs  are  in  contact  with  the  seat  pan  and  where  his 
head  is  in  contact  with  the  back  of  the  seat.  All  of  these  represent  sites 
of  potentially  serious  windblast  and  flail  injury.  Moreover,  Based  on  these 
results,  one  is  therefore  encouraged  to  pursue  the  mathematical  analysis  with 
a certain  degree  of  assurance  that  the  approach  is  faithfully  describing 
events  as  they  have  been  experimentally  observed  to  take  place.  With  this 
in  mind,  we  now  proceed  to  modify  the  theory  to  include  flow  separation 
effects . 

N 

SECTION  III 

VORTEX  MOTION  IN  A REGION  OF  SEPARATION 

In  aerodynamic  flow  theory,  two  simplifications  have  been  conveniently 
and  commonly  employed  to  Investigate  analytically  separated  flows  over  arbi- 
trary three  dimensional  bodies  of  revolution  (Marshall  and  Deffenbaugh,  1975). 
The  first  asserts  that  there  is  a direct  analogy  between  three-dimensional 
steady  flow  and  two-dimensional  unsteady  flow  (Allen  and  Perkins,  1951)  — 
such  that  a three-dimensional  steady  separated  flow  problem  can  be  analyzed 
as  a two-dimensional  unsteady  separated  flow  problem.  The  second  is  based 
upon  the  observation  (verified  experimentally)  that  a two-dimensional  unsteady 
wake  can  be  described  by  a distribution  of  inviscid  vortices,  originating 
from  the  separation  of  shear  layers  and  superimposed  on  the  unseparated  poten- 
tial flow  solution  (Mello,  1959,  Sarpkaya,  1968  and  Marshall  and  Deffenbaugh, 
1975).  These  vortices  are  modified  by  diffusion.  It  is  the  second  of  these 
two  assumptions  that  we  intend  to  pursue  here,  since  arguments  have  already 
been  presented  for  considering  flow  over  the  forearm  in  contact  with  a 
restraining  surface  to  be  two-dimensional  (Schneck,  1976). 
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Thus,  to  examine  the  separated  flow  patterns  around  the  geometric 
configuration  illustrated  in  Figure  2,  consider  the  situation  depicted  in 
Figure  6.  A vortex  pair  has  been  placed  on  the  lee  or  downstream  side  of 
the  two  cylinders  as  shown.  Neither  the  location,  yQ  + zq,  nor  the  strength, 
k,  of  these  two  vortices  can  be  specified  at  this  time.  In  fact,  the  location 
of  the  vortex  pair  can  not  be  determined  uniquely  from  potential  flow  theory 
alone.  One  must  either  consider  the  variation  of  pressure  across  the  boundary 
layer,  as  described  by  Sarpkaya  (1968)  and  Bar-Lev  and  Yang  (1975),  or,  one 
must  perform  a series  of  experiments  to  determine  empirical  values  for  y o 
and  zq.  The  intent  of  this  work  has  been,  and  is  to  exploit  both  of  these 
techniques  and  we  have  begun  with  the  latter. 

Depicted  in  Figure  7 is  a Scott-Engineering  Sciences  Model  9093  Hydro- 
dynamic  System  for  Flow  Visualization  Studies.  The  system  is  designed  to 
utilize  the  hydrogen-bubble  method  for  flow  visualization,  but  we  have  per- 
formed experiments  thus  far  with  a Pliolite  brand  of  solid  particles.  The 
Model  9093  provides  a self-contained  light  source  and  a simple,  but  effective, 
light  guide  which  enables  the  user  to  illuminate  almost  any  part  of  the  flow 
channel . 

In  addition  to  providing  means  for  observing  many  features  of  incompress- 
ible flow.  Model  9093  can  also  be  used  to  model  mathematically  a wide  range  of 
compressible  fluid  flow  behavior.  To  do  this,  the  well  known  analogy  between 
compressible  and  open  channel  flows  is  employed.  This  analogy  has  inevitable 
limitations  which  we  are  in  the  process  of  evaluating,  but,  nevertheless, 
there  are  experiments,  both  quantitative  and  qualitative,  which  are  worth 
performing.  Shown  in  Figure  8 are  some  early  results  which  we  have  obtained 
using  the  Model  9093  water  table,  into  which  was  placed  a doubel-cylinder 
configuration  such  as  that  depicted  in  Figure  2.  Note  the  vortices  forming 

on  the  downstream  side  of  this  geometrical  configuration  as  the  flow  moves 
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from  left  to  right.  Such  experimental  visualization  techniques  can  readily 
yield  information  relevant  to  the  further  development  of  the  mathematical 
theory  and  we  hope  to  be  able  to  continue  the  work  with  future  AFOSR  support. 
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